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Radio-loud AGN typically manifest powerful relativistic jets extending up to millions of light 
years and often showing superluminal motions organised in a complex kinematic pattern. A num¬ 
ber of physical models are still competing to explain the jet structure and kinematics revealed 
by radio images using the VLBI technique. Robust measurements of longitudinal and transverse 
velocity field in the jets would provide crucial information for these models. This is a difficult 
task, particularly for transversely resolved jets in objects like 3C273 and M87. To address this 
task, we have developed a new technique for identifying signihcant structural patterns (SSP) of 
smooth, transversely resolved flows and obtaining a velocity field from cross-correlation of these 
regions in multi-epoch observations. Detection of individual SSP is performed using the wavelet 
decomposition and multiscale segmentation of the observed structure. The cross-correlation algo¬ 
rithm combines structural information on different scales of the wavelet decomposition, providing 
a robust and reliable identification of related SSP in multi -epoch images. The algorithm enables 
recovering structural evolution on scales down to 0.25 full width at half maximum (FWHM) of the 
image point spread function (PSF). We present here examples of applying this algorithm to obtain 
the first detailed transverse velocity fields and to study the kinematic evolution in the parsec-scale 
jets in 3C 273 and M 87. 
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1. Introduction 

The steady improvements of dynamic range of astronomical images and increasing complex¬ 
ity and detail of astrophysical modeling bring a higher demand on automatic (or unsupervised) 
methods for characterization and analysis of structural patterns in astronomical images. 

A number of approaches developed in the fields of computer vision and remote sensing for 
tracking structural changes [c/, 1] typically require oversampling in the temporal domain. This 
renders them difficult to be used in astronomical applications that requires tracking changes moni¬ 
tored with sparse sampling and dealing with partially transparent, optically thin structures. 

Presently, structural decomposition of astronomical images normally involves simplified su¬ 
pervised fechniques based on idenfificafion of specific feafures of fhe sfrucfure [2] or fiffing fhe 
observed sfrucfure wifh a sef of predefined femplafes {e.g., fwo-dimensional Gaussian feafures). 

A more robusf approach fo automatize idenfificafion and fracking of sfrucfural pafferns in 
asfronomical images can be provided by a generic mulfiscale mefhod such as wavelef deconvolution 
or wavelef decomposifion [cf., 4]. In [5], we explored fhis wavelef approach and presenf a wavelef- 
based image segmenfafion and evaluafion (WISE) mefhod for sfrucfure decomposifion and fracking 
in asfronomical images. In fhis confribufion, we presenf an overview of fhis new mefhod along wifh 
resulfs obfain on WISE analysis of parsec-scale radio jefs in 3C273 and 3C 120 compared wifh 
resulfs of convenfional sfrucfure analysis previously applied fo fhese dafa. Einally, preliminary 
resulfs of fhe velocify field analysis of fhe jef in M 87 is presenfed. 

2. Wavelet-based Image Segmentation and Evaluation 

To characferize fhe sfrucfure and sfrucfural evolufion of an asfronomical objecf, fhe imaged 
objecf sfrucfure needs fo be decomposed info a sef of significanf sfrucfural pafferns (SSP) fhaf can 
be successfully fracked across a sequence of images. The mulfiscale decomposifion provided by fhe 
wavelef fransform [6] makes wavelefs exceptionally well-suifed fo perform such a decomposifion, 
yielding an accurafe assessmenf of fhe noise variafion across fhe image and warranfing a robusf 
represenfafion of fhe characferisfic sfrucfural pafferns of fhe image. To adapf beffer fo fhe needs 
of asfronomical imaging, we have exfended fhe mulfiscale approach fo objecf defecfion, similarly 
fo fhe mefhodology developed for fhe mulfiscale vision model [MVM; 4]. The mefhod employs 
segmenfed wavelef decomposifion (SWD) of individual images info arbifrary fwo- dimensional 
SSP (or image regions) and subsequenf mulfiscale cross-correlation (MCC) of fhe resulfing sefs of 
SSP. 

2.1 Segmented Wavelet Decomposition 

The segmenfed wavelef decomposifion (SWD) comprises fhe following sfeps fo describe an 
image sfrucfure by a sef of sfafisfically significanf pafferns: 

1. A wavelef fransform is performed on an image I by decomposing fhe image info a sef of J 
sub-bands (scales), wj, and esfimafing fhe residual image noise (variable across fhe image). 

2. Af each sub-band, sfafisfically significanf wavelef coefficienfs are exfracfed from fhe decom¬ 
posifion by fhresholding fhem againsf fhe image noise. 
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3. The significant coefficients are examined for local maxima, and a subset of the local maxima 
satisfying detection criteria is identified, defining the locations of SSP in the image. 

4. Two-dimensional boundaries of the SSP are defined by the watershed segmentation using 
the feature locations as initial markers. 

The decomposition provides a structure representation that is sensitive to compact and marginally 
resolved features as well as to structural patterns much larger than the FWHM of the instrumental 
PSF in the image. 

2.2 Multiscale Cross-Correlation 

To detect structural differences between two images of an astronomical object made at epochs 
ti and t 2 , one needs to find an optimal set of displacements of the original SSP (described by the 
groups of SSP Sjj,j = 1, ...,7) that would match the SSP in the second image (described by Sj^ 2 , 
j = 1,...,7). Cross-correlating Sj^i and 5'y2 is a well-suited tool for this purpose. There are two 
specific issues that need to be addressed, however, to ensure that the cross-correlation analysis 
is reliable. First, a viable rule needs to be introduced to identify the relevant image area across 
which the cross-correlation is to be applied. Second, the probability of false matching needs to be 
minimized for features with sizes smaller than the typical displacement between the two epochs. 

These two requirements can be met by multiscale cross-correlation (MCC), which combines 
the structural and positional information contained in Sj at all scales of the wavelet decomposition: 

1. The largest scale 7 of a wavelet decomposition is chosen such that the largest expected dis¬ 
placement is smaller than 2“^. 

2. Displacements of SSP features are determined at the largest scale J. For this calculation, all 

are set to zero, and ^ = 5j^a is calculated for each SSP. 

3. At each subsequent scale j {j < J), Aj^J is determined first by adopting the displacement 
Aj^a measured at the j 1 scale for the SSP in which the given y-scale region sj^a falls. Then 
the total displacement for this SSP is given by Ay ^ = Aj^J -|- dj^a- 

In this algorithm, the only quantity that needs to be calculated at each scale is the relative displace¬ 
ment 5j^a- This quantity can be determined reliably from the cross-correlation. 

3. Applications to astronomical images 

We have tested the performance of WISE on astronomical images by applying it to several 
image sequences obtained as part of the MOJAVE long-term monitoring program of extragalactic 
jets with very long baseline interferometry (VEBI) observations [7, and references therein] and in 
the M 87 “movie” project [8]. 

The analysis of the jet in 3C 273 was performed on 69 images, with the observations covering 
the time range from 1996 to 2010 and providing, on average, one observation every three months. 
The core separations of individual SSP obtained from WISE decomposition are compared in Eig. 1) 
with the results from the MOJAVE kinematic analysis based on the Gaussian model fitting of the jet 
structure. Comparison indicates that WISE detects consistently nearly all the components identified 
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Figure 1: Core separation plot of the most prominent features in the jet of 3C 273 (taken from [5]). The 
model-ht based MOJAVE results (dashed lines) are compared with the WISE results (solid lines). 


by the MOJAVE model fitting analysis, with a very good agreement on their positional locations 
and separation speeds. 

Similar analysis was performed on 87 observations of the jet in 3C 120 between 1996 and 
2010. A total of 30 moving SSP was detected. The resulting core separations of the SSP plotted 
in Fig. 2 generally agree very well with the separations of the jet components identified in the 
MOJAVE Gaussian model fit analysis. 

The jet in M 87 is probably one of the most investigated. Its proximity (D = 16.7 Mpc; [9]) 
combined with a large mass of the central black hole (Mbh — 6.1 x IO^Mq; [10]) make it one of the 
primary source to probe the jet formation and acceleration at its base. We analyzed eleven epochs 
of the M 87 jet observed with the Very Earge Baseline Array (VEBA) between 27 January 2007 
and 26 August 2007 with an average of 21 days interval [8]. The velocity field obtained from WISE 
analysis is plotted in Fig.3. It reveals complex kinematic with subluminal and superluminal motion 
measured in three regions of the jet: a southern and northern sheath and a central spine. 

A Stacked Normalized Cross Correlation analysis combining all SSPs obtained at all eleven 
epochs was performed. Result indicates the presence of 2 main speed component in all three 
regions of the jet. A slow mildly relativistic speed (~ 0.5 c), and a faster relativistic speed (~ 2.25 
c). This could be explained by the presence of a disk wind around the jet from which only a small 
fraction of it at the interface could be observable thi'ough synchrotron radiation [11]. 

Additionally we found that the average velocity measured in the southern sheath was consis¬ 
tently faster than the average velocity measured in the northern sheath, suggesting jet or pattern 
rotation. 
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Figure 2: Core-separation plot of the features identified in the jet of 3C 120 (taken from [5]). The model-fit 
based MOJAVE results (dashed lines) are compared with the WISE results (solid lines). 
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Figure 3: Velocity field in the jet of M 87. Three main regions are detected: a southern (blue) and northern 
(green) sheath and a central spine (red). Arrows represents displacement of individal SSPs between two 
subsequent epochs. 
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4. Conclusions 

The WISE method we presented here offers an effeetive and objeetive way to elassify strue- 
tural patterns in images of astronomieal objeets and traek their evolution traeed by multiple ob¬ 
servations of the same objeet. The method eombines automatie segmented wavelet deeomposition 
with a multiseale eross-eorrelation algorithm, whieh enables reliable identifieation and traeking of 
statistieally signifieant struetural patterns. 

We performed the analysis of several transversely resolved AGN jets observed as part of the 
MOJAVE projeet using this method. Eor all this sourees, global kinematie ehanges obtained using 
our SWD teehnique were in exeellent agreement with the one obtained by the MOJAVE team using 
model-fitting teehnique, demonstrating the reliability of our method for the reeonstruetion of the 
veloeity field in transversely resolved flows. 

Preliminary results of the applieation of WISE on the jets of M 87 suggests a layered jet with a 
fast spine surrounded by a subrelativistie disk wind. Ongoing analysis will eoneentrate on attempt¬ 
ing to model veloeity in eombination with the eollimation profile, provide an hint on the potential 
jet rotation, and investigate the veloeities obtained in the eounter jet side. 
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